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FORMULAS  FOR  SOURCE,  DOUBLET,  AND  VORTEX  DISTRIBUTIONS 

IN  SUPERSONIC  WING  THEORY 

By  Harvard  Lomax,  Max.  A.  Heaslet, 
and  Franklyn  B  .  Fuller 

SUMMARY 


The  formulas  of  supersonic  wing  theory  for  source,  doiiblet,  and 
vortex  distributions  are  reviewed  and  a  systematic  presentation  is 
provided  which  relates  these  distributions  to  the  pressure  and  to  the 
vertical  induced  velocity  in  the  plane  of  the  wtug.  It  is  shown  that 
care  must  be  used  in  treating  the  singularities  involved  in  the  analysis 
and  that  the  order  of  integration  is  not  always  reversible.  Further,  it 
is  shown  that  the  use  of  the  complex  variable  can  often  facilitate  the 
calculation  of  the  integrals  involved.  Certain  special  applications  are 
included  to  illustrate  the  concepts  presented. 


INTRODUCTION 


One  of  the  most  fundamental  approaches  to  the  analytical  investiga¬ 
tion  of  linearized  wing  theory,  throughout  the  subsonic  and  supersonic 
Mach  nuinber  range,  stems  from  the  use  of  certain  elementary  mathematical 
expressions  which  are  identified  physically  with  sources,  doublets,  and 
vortices  in  the  fluid  medium.  By  means  of  these  expressions,  boundary- 
value  problems  involving  wings  with  thickness,  camber,  and  angle  of 
attack  can  be  solved.  These  problems  are  divided  into  two  categories: 
one,  involving  symmetrical  bodies  with  thickness  and  no  lift,  is  analyzed 
by  means  of  source  distributions;  and  the  other,  involving  lifting  plates 
without  thickness,  is  analyzed  by  means  of  doublet  and  vortex  distribu¬ 
tions. 

) 

All  these  distributions  require  the  treatment  of  singuleLrities  in 
the  mathematical  analysis.  Thus,  for  subsonic  Mach  numbers,  the  concept 
of  Cauchy*  s  principal  part  plays  an  inqportant  role  in  the  calculation  of 
integrals  arising  in  the  development  of  lifting— line  and  two-dimensional 
section  problems.  In  supersonic  wing  theory,  the  Cauchy  principal  part 
is  again  used  in  the  treatment,  for  example,,  of  conical— flow  problems  as 
in  reference  1,  but,  because  of  the  Mach  lines  and  cones  appearing  in 
the  physical  flow  and  in  the  hyperbolic  geometry  of  the  differential 
equation,  other  techniques  in  handling  improper  integrals  are  needed. 
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The  integrals  in  supersonic  wing  theory  thus  require^  in  general^  more 
carefial  attention  to  the  discontinuities  in  the  integrand  and^  as  an 
illustration^  indiscriminate  use  of  such  standard  devices  as  inversion 
of  the  order  of  integration  may  lead  to  incorrect  results. 

When  problems  of  the  first  kind  are  involved,  that  is,  when  pre¬ 
scribed  distributions  are  to  be  integrated  (as  for  the  problem  of  finding 
the  pressure  on  a  wing  with  symmetrical  thickness),  a  guide  to  the  proper 
method  of  calculation  is  often  furnished  by  physical  intuition.  However, 
when  problems  of  the  second  kind  arise,  that  is,  problems^  the  solution  of 
which  depends  upon  the  inversion  of  an  integral  equation  (as  the  flat 
plate  of  arbitrary  plan  form),  the  mathematical  methods  Eire  more  abstract. 

The  purpose  of  the  present  report  is:  first,  to  review  the  formulas 
of  linearized  wing  theory  in  which  source,  doiiblet,  and  elementary- 
vortex  distributions  are  introduced  and  to  relate  these  distributions  to 
the  pressure  and  to  the  vertical  induced  velocity  in  the  plane  of  the 
wing;  second,  to  show  that  the  use  of  the  complex  variable  can  often 
facilitate  the  calculation  of  the  integrals  involved;  and  finally,  to 
present  certain  special  applications  which  will  illustrate  the  basic 
concepts. 


LIST  OF  IMPORTANT  SYMBOLS 


Mo 

q 

q 

r„ 


pressure  coefficient 
free— stream  Mach  number 

loading  coefficient  (pressure  on  the  lower  surface  minus  pressure 
on  the  upper  surface  divided  by  free— stream  dynamic  pressure) 


free— stream  dynamic  pressure 


1  Po 


hyperbolic  distance  between  points  x,y, z  and  xi,yi,0; 
r^  =J  (x-xi)^ 


Yq  velocity  of  the  free  stream 
x,y,z  Cartesian  coordinates 
r,s,z  characteristic  coordinates 
|,T|,z  oblique  coordinates 

mi  cotangent  of  the  angle  between  the  and  x  axes 
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3 


ms 

^2 

U 

V 

p 

X 

9 

A 


cotangent  of  the  angle  between  the  |  and  x  axes 
J 1+mi® 
y l+ms^ 

perturbation  velocity  in  x  direction 
perturbation  velocity  in  2  direction 

*  C 

slope  of  stream  surface  ('w/Yq) 
perturbation  velocity  potential 

Jump  in  value  of  the  quantity  considered  across  the  z  =  0  plane 


Subscript 

u  value  of  a  quantity  on  the  upper  surface  of  a  wing  (z  =  0  plane) 


\ 

GENERAL  THEORY 


It  is  well  known  that  the  analysis  of  thin  wings  at  supersonic 
speeds  and  at  small  angles  of  attack  can  be  expressed  in  mathematical 
terms  as  a  boundary— value  problem  for  the  wave  equation.  If  Q  repre¬ 
sents  a  velocity  potential  or  any  one  of  the  velocity  components  them¬ 
selves,  this  equation  can  be  written 


p  ^xx  ^yy  ^  ^zz  ^  (l) 

where  the  z  =  0  plane  is  the  plane  of  the  wing,  the  free— stream  veloc— 
Ity  Yq  is  directed  along  the  x  axis  and  =  Mo^-"l>  Mq  being  the 
free— stream  Mach  number. 

Of  the  many  ways  of  solving  the  boundary— value  problems  associated 
with  the  wave  equation,  the  most  convenient  for  the  present  purpose  is 
the  Yolterra  solution.  In  reference  1,  a  discussiozi  was  given  of  the 
application  to  aerodynamic  problems  of  Yolterra*  s^^thod.  Thus,  the 

solution  for  in  terms  of  A  ft  and  ^he  jump  in  ^  and  its 

gradient  in  crossing  the  z  =:  0  plane  canl>e  written,  if 


k 
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rc  =J  (x-xi)^  (y-yi)^  -P^z^ 


in  the  form 


i_  A 

25t  Sx 


arc  cosh 


PvCy-yi)^  + 


dxi  dyi  + 


2k  Sx 


z(x-xi)  dxi  dyi 
[(y-yi)2  +  z^]  r, 


where  the  region  t  is  that  portion  of  the  z  =  0  plane  lying  within 
the  forecone  from  the  point  x,y,z  (the  forebranch  of  the  hyperbola 

rc^=0). 

Equation  (3)  represents  a  general  solution  to  the  yave  equation  and 
has  yet  to  be  put  in  a  form  which  represents  directly  a  solution  to  a 
problem  arising  in  the  study  of  wings.  When  equation  (3)  is  so  adapted 
in  the  next  section^  it  will  represent  the  velocity  potential  due  to  a 
distribution  of  sources^  doublets^  and  elementary  horseshoe  vortices, 
the  strength  of  which  are  given  in  terms  of  the  wing  shape  and  loading. 

The  following  study  of  the  adaptation  of  equation  (3)  to  the  partic¬ 
ular  boundary— value  problems  of  wing  theory  requires  the  introduction  of 

two  axial  systems  other  than 
the  x,y,z  Cartesian  coordi¬ 
nates  already  defined.  First, 
the  5^T|,z  coordinate  system 
is  defined  so  that  z  is  normal 
to  the  plane  of  the  wing  while 
y-f^  bTi  are  both  normal  to  z 

lie  in  the  plane  of  the 
y  wing)  and  make  arbitrary  angles 

x=  with  the  x,y  system  (see  sketch) 


^  M-2  =  V  I+212' 


x=miy 


The  equations  which  relate  the 
|,Ti,z  to  the  x,y,z  system 
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7 


z 


^  42(x-may) 

mi+m2 

=  JL_  J_ 

T,  =  4i(z+mgy) 

mi+mg 

=  Z 

z  =  z 

(5) 


and  the  Jacohidn,  the  relation  between  the  differential  areas  in  the  two 
systems,  is 


mi+me 

dz  dy  =  -  d|  dTi 


(6) 


Finally,  the  value  of  r^  is  transformed  by  the  equation 


di‘ 


^  2(n-^li)(g-gi)(miing+P^) 


*  ,.2 


(7) 


Second^  the  r^s^z  coordinate 
system  is  also  defined  so  that 
z  is  normal  to  the  plane  of 
the  wing  while  r^s  are  both 
normal  to  z  and  lie  along 
the  traces  of  the  Mach  cone 
emanating  from  the  origin 
(see'  sketch).  It  is  apparent 
that  the  r^s^z  coordinates 
are  a  special  case  of  the  |,t),z 
system  formed  when  mi=m2=3. 
(Notice  also  that  the  x,y,z 
coordinates  are  obtained  from 
^r\,z  when  mi=0  and  1112=00.) 
The  equations  which  relate 
r,s  and  z  to  z,y  and  z 
are 
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^=A(s+r) 

y.i(s-r)  =-|(=‘+Uy) 


(8) 


(9) 


(10) 


THE  'I'H'r™  FDKDAMEIOTAL  FOBMUIAS 


As  has  heen  indicated  already,  the  next  purpose  is  to  relate  equa¬ 
tion  (3)  to  the  three  fundamental  formulas  arising  in  wing  theory  which 
are  those  relating  the  velocity  potential  to  source,  doublet,  and  vortex 
distrlhutions.  These  distributions  can  be  expressed  in  terms  of  the 
discontinuity  in  either  cp  or  its  gradients  in  the  plane  of  the  wing. 

In  this  way,  the  expression  for  a  source  distribution  is  obtained  when 
the  perturbation  velocity  potential  9  is  an  even  function  with  respect 
to  the  z  =  0  plane  and  is  expressed  as  a  doxible  integral  involving 
A(^9/5z)  where  the  A  notation  denotes  the  jump  in  the  value  of 
Scp/Sz  in  crossing  this  plane;  that  is 

hz  -4=0+  >'z=0- 

Similarly,  a  distribution  of  doublets  is  obtained  when  9  is  an  odd 
function  with  respect  to  the  z  =  0  plane  and  is  expressed  as  a  double 
integral  involving  A  9.  Finally,  a  vortex  distribution  results  when  9 
is  an  odd  function  and  given  as  a  double  integral  involving  the  loading 
coefficient  Ap/q  which  is,  in  turn,  equal  to  2A9x/Yo*  these 
formulas  A(Wdz),  A9,  and  Ap/q  are,  respectively,  the  strengths 
per  imit  area  of  the  sources,  doublets,  and  elementary  horseshoe  vortices. 
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The  Source  Distribution 


The  velocity  potential  at  the  point  x.,y,z  of  a  unit  source  at  the 
point  xi,yi,0  is  given  by  the  equation  cp  =  —  l/Sitr^,  The  manner  in 
which  a  distribution  of  these  sources  affects  the  potential  will  now  be 
developed  from,  equation  ( 3)  •  If  >  in  equation  ( 3)  ^  is  set  equal  to 
the  velocity  potential,  and  the  potential  in  turn  is  assumed  to  be 
symmetrical  above  and  below  the  x,y  plane  (as  in  the  case  of  a  symmet¬ 
rical  airfoil  at  zero  lift),  then  equation  (3)  becomes 


fp  =  - 


1  b 


2n  dx 


A 


d(p 


arc  cosh 


X— Xl 


py(y-yi)^+z^ 


dxi  dyi 


Since  the  inverse  hyperbolic  term  vanishes  on  the  Mach  forecone,  the 
partial  derivative  can  be  carried  through  the  double  integral  sign  and 
there  results 


cp  (x,y,z)  =  -  i 


^c 


dxi 


dyi 


(11) 


vhere  is  the  vertical  induced  velocity  on  the  upper  side  of  the 

z  =  0  plane  and  A(Scp/Sz)  =  2wu.  by  reasons  of  symmetry.  Equation  (ll) 
is  the  familiar  equation  for  the  velocity  potential  due  to  a  distribution 
of  sources  in  the  xy  plane. 

In  the  coordinate  system  equation  (11)  becomes 

^  - - /  / - d^i  df]!  (12) 

where  the  area  t  is  transferred  to  the  plane,  and  r^,  is  given 

in  these  coordinates  by  equation  (T)« 


The  Doublet  Distribution 


The  velocity  potential  at  the  point  x,y,z  of  a  doublet  at  the 
point  xi,yi,0  is  given  by  the  equation  cp=  P^z/2jcrj,®.  The  effect  on 
the  potential  of  a  distribution  of  these  doublets  is  not  so  obvious  as 
it  was  in  the  case  of  the  sources  and  considerable  care  must  be  shown 
in  the  development. 
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If  the  vertical  induced  velocity  is  considered  always  to  be  equal 
above  and  below  the  xy  plane,  but  the  potential  to  be  discontinuous 
across  it,  then  equation  ( 3)  for  the  velocity  potential  becomes 

[(y-yi)2+z2]  r^ 


In  this  case  the  integrand  does  not  vanish  at  the  Mach  cOne  and  the 
partial  derivative  cannot  be  moved  directly  through  the  double  integral 
sign.  Writiiig  in  the  limits  of  integration  so  that  the  first  integration^ 
is  made  with  respect  to  yi, 

T  nx-^z  pYa  z(x-xi)  A  cp  (xi,yi) 

cp  =  ^  |-  /  dxi  /  dyi  - ,  . . . -  . . 

oxJ«„  Jyi  [(y-yi)^+z^]  y(x-xi)^-P^(y-yi)^-P^z^ 

(13a) 

vhere 

Yi  =  y  -  (x-xi)^  -P^z^ 

and 

Ta  =  y  +  y(x-xi)^  -P^z^ 

If  In  the  yi  Integral  xi  is  replaced  by  the  value  x-Pz,  the 
result  is  indetermi3aate .  Such  an  indeterminate  form  can  be  evaluated, 
however,  by  excluding  the  limit  x-Pz  from  the  area  of  integration. 

Hence  consider  the  integral 


_ z(x-xi)  A(p(xi,yi) _ 

[(y-yi)^+z^]  y (x-xi)®-P^(y-yi)^-^®z^ 


^Since  in  what  follows  the  order  of  Integration  is  important,  the  nota¬ 
tion  will  be  adopted  that  /dy  /clx  f(x,y)  =  /  [/f(x,y)dx  Jdy;  that  is, 
the  integration  is  made  first  with  respect  to  x.  When  the  notation 
ff  f(x,y)dx  dy  is  used,  the  order  of  integration  is  Immaterial. 


NACA  TN  2252 


By  formal  manipulation 


=  lii 


py+€  yz^+e^  A  cpCx-^yz^+e®,  yi)dyi 

T-e  [(y-yi)^+z®]  yeMy-yi)® 


lim  1  rx-p  y  z^+e^  a  pTa 


_ 2(x-xi)  A9  (xi,yi) _ 

[  (y-yi)^+z^]  y  (x-Xi)®-P®(y-yi)®-^®z^  • 


(13b) 


Application  of  the  mean  value  theorem  to  the  first  Integral  in  equa¬ 
tion  (l3b)  yields 


lim  —  V  z^+e^  A  <p(x-P  V  z^+e^,  y+06) 

€— >0  2jt 


y-e  [  (y-yi)^+z^]y €My-yi)^ 


which  becomes 


^0  2  ‘‘^  ^(x-pyz^+e^,  y+0€)  = 


I  Acp(x-0z,y) 


The  second  Integral  in  equation  (l3b)  is  simplified  by  Introducing 
the  notation  of  the  finite  part.  Define  the  symbol  -f  by  the  eqiiation 


f(y)<3y  ^  _  1  ^  f(y)dy  ^  ^ 

o  Sz  V  ei-y  2  hx'^o  (a— y)®/^  Sx 


f(y)dy 

Va-y 


When  applied  to  a  single  integral  this  definition  is  consistent  with  the 
usual  ones  for  the  finite  part  given  as 
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f(y)dy 
o  (a-yf/s 


f(y)-f(a) 
'^o  .(a-yf/2 


f(a) 

V“i“ 


f*(y)<iy 
V  a-y 


+ 


\ 


However,  when  applied  to  double  integrals  an  inconsistency  with  regard  to 
the  order  of  integration  between  the  two  symbols  |  and  ^  arises. 
Hadamard  (reference  2)  and  Eoblnson  (reference  3)  both  use  the  convention 

that  the  order  of  integration  in  the  operation  [  //  f(z,y)dy  dx  is 
reversible;  that  is, 


/  dy  /  dz  f(x,y) 


/  dx  /  dy  f(x,y) 


Such  a  convention  excludes  from  the  area  of  integration  all  singularities 
over  which  the  order  of  integration  is  not  reversible.  These  singular 
regions  are  then  treated  separately.  This  convention  has  the  disadvantage 
that,  in  constructing  a  series  of  integrals,  the  value  of  a  given  inte¬ 
gral  is  not  independent  of  succeeding  integrals. 

The  operator  -f  avoids  the  above  difficulty.  The  value  of  an  inte¬ 
gral  defined  by  f(y)dy  is  independent  of  succeeding  operations.^  At 
the  same  time,  however,  the  order  of  integration  of  operations  involving 
the  sign  cannot  be  reversed.  Hence 

f  dy  ^dx  f(x,y)  ^  /dx  -/-dy  f(x,y) 


^For  exeunple,  according  to  reference  2 


ri  dll 

'o  (ni)°^=v^ 


but  according  to  the  same  reference 


However, 


=  -  2Jt 


^ _ r  de  _ 
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Eetum  now  to  the  second  integral  in  equation  (l3b).  Applying  the 
aefinition  already  mentioned  for  f-,  it  is  found  for  this  type  of  inte¬ 
gral  that 


^  nA(x)-K}(x)  ay  ^ 

A(x)-G(x)  yG^-(A-y)^  J 

By  means  of  the  last  formula^  the  equation  for  9  hecomes 

cp  =  i  Acp(x-pz,y)- 1^  J  dxi  ^  dyi  (l4a) 

where  r^  is  given  by  equation  (2). 

Notice  that  if  the  integration  had  been  made  first  with  respect  to 
xi  the  result  would  be 


1  ^ 


f(x,y) 


Sz  LyGMA-y)2  J 


dy 


_ z(x-xi)A  9(^i^yi) _ 

[ ( y-yi ) ] y ( Z-Xi ) ( y-yi ) ^-P^z^ 

which  reduces  immediately  to 

(I4t) 

J  J  q-  Tq 

/ 

Equations  (l4a)  and  (l4b)  illustrate  the  vital  importance  of  the 
order  of  integration.  In  fact  by  subtraction  the  result  can  be  derived 
that 


dxi  7  dyi 


A  cp 
,  3/2 


i 


dyi  4-  dxi 


A  9 

TeTi 


zp‘ 


A9  (x-Pz,y) 


If  the  coordinates  are  used,  equation  (13)  becomes  in  the 

notation  presented  in  equation  (5) 
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_i_  A  +  1. 

2n  \|i2  Hi  hij 


z  (n-'li)  —  +  (l-li)  —  AcpdT^idlj^ 
L  Hi  ^^2  J 


(ti-^i)  (l-li)  ^  +z2[r, 

^1  H2  J  J 


The  limits  of  integration  depend 
now  upon  the  position  of  the 
axes  I  and  with  respect  to 
the  Mach  lines  in  the  x,j  plane. 
The  area  t  is  still  hounded  by 
the  curve  rg^=0  and  infinity. 

The  asymptotes  for  the  curve 
Tq^=0  are  given  by  the  two  equa¬ 
tions 


ma  <  P 


(16) 


The  sketch  shows  how  the 
area  t  in  the  xy  plane  trans¬ 
fers  to  the  |,Ti  plane.  In  the 
case  in  which  both  mi  and  m2 
are  less  than  P  (the  case  for 
which  the  sketch  was  drawn) ,  the 


n^i  = 


=  -— (l-ii) 


asymp'bo'bes  are  sl^raiglrt  lines  having  positive  slopes  and  the  limits  of 
integration  are  always  from  one  side  of  the  cone  to  the  other  and  from 
minus  infinity  to  some  TnftT-iTmmi  value*  Consider  the  case  in  which  the 
integration  is  made  first  with  respect  to  t|i*  Then  defining 


there  results 
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which,  reduces  to 


where 


I  A  fp  (6a>na) 


Tla 


z^2  V 
ini+iii2 

zdi(iiiime+3^) 
(mi+nte)  J 


(20) 


} 


second  term  in  equation  (19)  simplifies  when  the  finite— part  notation 
is  introduced  so  that  finally 


cp  =  |A(p(la^na) 


zP^Ciiii+me) 

2Jt|li|l2 


.  Acp(Ei^Tli) 
-iTli  - j:j3  — 


(21) 


If  the  Integration  had  been  made  first  with  respect  to  the  results 

would  have  been 


cp  =  I A  cp(  - 


zp^(mi+mg2 

27t}J.iU2  J  Jr 


1  .  r  ^ 


(22) 


where 


%  =  n 


ziii  y 

mi+ms 


(23) 


(mi +1112) 


The  geometric  interpretation  of  the  points 
as  follows:  They  are  the  points  at  vhich  the  forecone  in  the  |  yr\ 
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plane  (given  "by  one  ‘brancli  of  the  curve  attains  maximum  values 

for  r\  and  respectively,  (see  sketch). 


As  the  g  and  r\  axes  approach 
the  Mach  lines  in  the  x,y  plane,  that 
is,  as  m2  and  mi  approach  p,  the 
residue  terms  in  equations  (2l)  and  (22 
approach  (1/2)A  cp(|a^"®°) 

(1/2)A  cp(— oo,T]i^),  respectively,  which 
represent  the  jump  in  potential  infi¬ 
nitely  far  distant  from  the  point  P 
(and  hence  may  he  taken  as  zero).  Thus 
when  the  axes  lie  along  the  Mach 

lines,  thereby  becoming  the  r,s  axes 
of  equations  (3)>  the  equations  for  cp 
are  without  the  residue  terms  and  the 
order  of  integration  is  immaterial. 

When  mi  and  m2  are  greater  than  P 
the  same  is  true  (i.e.,  the  terms 
(l/2)A<p(|a,Tia)  and  (i/2)A 
are  missing  from  equations  (2l)  and 
(22),  respectively),  so  that  the  effect 
of  a  distribution  of  doublets  on  the 
velocity  potential  can  be  summarized  as 
being 


ce  of 
econe 
plane 


for  0<mi<p,  0<m2< 


'P  =  iA<p(S,„3)-2p^ 

2  “  2nM,i)i2 


<i§i  /  dill 


(24a) 


for  0<iii2<P,  0<mi< 


tP  =  i  Acp(ib^%)  - 


zP^(mi+ii)a) 

2jtM-i|J2 


(24b) 


for  P  <  mi  <  00,  P  <  ni£  < 


zp2(mi+iD2) 


Aq)(  ,Tii) 


dTiidl 


(24c) 


where  la^’la^  Ib^  given  by  equations  (20)  and  (23). 
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There  exists  the  interesting  corollary  obtained  by  subtracting  equa¬ 
tion  (24a)  from  (24b);  namely,  that  the  difference  between  an  integration 
of  supersonic  doublets  made  first  in  one  order  and  then  in  the  reverse  is 
equal  to  the  difference  in  the  magnitude  of  the  distribution  at  two  points 
in  the  plane. 


The  Vortex  Distribution 


The  velocity  potential  at  the  point  xi,yi,z  of  an  elementary 
horseshoe  vortex  dt  the  point  xi,yi,0  is  given  by  the  equation 
cp  =  -z(x-xi)/2Jt  [(y-yi)^+z^]rg.  If  0  in  equation  (3)  is  taken  to  be 
the  induced  velocity  u  in  the  direction  of  the  free  stream,  and,  as  in 
the  case  involving  doiiblets,  the  flow  field  is  considered  to  contain  no 
bodies  with  thickness  (so  A(^u/Sz)  is  everywhere  zero),  then  eqiia- 
tion  (3)  reads 


_1_  r  rz(x— xi)Au  dxidyi 

2it  J  [(y-yi)^+z^]rc 


Now,  since  by  definition 


the  potential  for  a  distribution  of  vortices  can  be  written 


(25) 


« 


1_  r  r  z(x-xi)Au(xi,yi)dxidyi 

2it  JvJ  [(y-yi)^  +z^] 


In  terms  of  the  |,ti,z  coordinates,  equation  (26)  transforms  by 
simple  substitution  to  the  equality 


cp  = 


(mi+m2)z 

2jtpi|j.2 


(.-ix)  ^  +  (Mx)  S 


dll  dqi 


(27) 
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APPLICATION  OF  FiraDAMENTAL  FOEMCJLAS  TO  THE  PLANE  OF  THE  WING 


Pressure  in  Terms  of  Slope  for  Symmetrical  Bodies 


Since  in  linearized  theory  the  equations  for  pressure  coefficient 
and  surface  slope  are 


the  equation  for  the  soxarce  distribution,  equation  (ll),  can  be  rewritten 
as 


2  h  r  rWxi,yi)dxidyi 


where  is  the  value  of  the  slope  on  the  upper  s-urface  of  a  symmetri¬ 

cal  body.  In  terms  of  the  |,tj,z  system  this  becomes  (compare  the 
transition  from  equation  (13)  to  equation  (15)) 


P  \  M-2  Sq  Hi 


^u(  tl^Tji)d|i  drii 


In  carrying  the  partial  derivative  through  the  integrals  there  results 
(as  in  equa,tlon  (19))  for  0  <  mi  <  p  the  equation 


Cp  =  lim 


/  zH2  //P^-mi^  \ 

>Lo*+Li»  I  ■  mi-f-me — ) 


0  «l^i^Lo'-Li« 


/  \  mi  /  V  mg 

2(mi+ms,)  p  p  (n^i)  ^7  +  (l-l^)  — 


and  this  reduces  to  the  expression 
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for  0  <]ni  <p,  ,0  <iii2  <00 


^u( 


2(mi+iii2) 


r  ^ 


'U 


(31a) 


Similarly,  for  0<m2<P,  0<mi<oo 


^u(|-b>Tl-j3)  - 


2(mi+mg) 


1^1 


/ 


(Sl-b) 


and  for  P<mi<oo,  P<iii2<eo 


'P 


2(mi+m2) 

«MiP2 


I L 


.  .  mi  ,  .  nfe 

(Ti-rii)  —  +  (l-l^)  — 


rc' 


<111  (31c) 


where  given  by  equations  (20)  and.  (23) 

As  special  cases  of  the  results  given  by  equations  (31),  consider  ^ 
and  T]  to  represent  the  x,y  and  the  r,s  axes,  respectively.  In  the 
former  case,  x  replaces  ^  as  ms  — ^  00  and  y  replaces  q  as 
mi  ->  Oj  hence 


Cp  =  I  Mx-Pz,y)  -  I J' dxi^dyi 


(x-xi)Xu(xi,yi) 


(32a) 


and 


(x-xi)^u(xi,yi) 


(32b) 


In  the  latter  case  r  replaces  |  as  and  s  replaces  t]  as 

Dll  — ^  P;  hence 
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In  each  of  the  equations  (3l)  throng  (33)#  z  can  he  set  equal  to 
zero  without  affecting  the  validity  of  the  equalities.  When  z  vanishes# 
the  following  identities  hold 


With  these  equations  the  pressure  coefficient  on  the  surface  of  a  syai- 
metrical  nonlifting  wing  can  be  determined  if  the  surface  slope  is  given. 
The  special  cases  of  equations  (3l)  in.  the  plane  z=0  and  in  the  x#y 
and  r#s  coordinate  systems  are  given  in  the  summary  of  this  section 
(see  equations  (4?)  and  (50)). 

Vertical  Induced  Velocity  in  Terms  of  the  Jump  in  Potential 

It  is  proposed  next  to  find  the  vertical  induced  velocity  in  the 
z=0  plane  as  a  function  of  the  jump  in  potential  across  that  plane. 

Consider  equations  (24)  for  the  doublet  distribution  and  take  the 
partial  derivative  with  respect  to  z  of  both  sides;  then  find  the  limit 
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of  the  resulting  expression  as  z  goes  to  zero.  If  equation  (24a)  is 
used, for  example,  there  results  for  the  first  term 


lim 
z  0 


a  1 
2 


Acp  (1 


2  Sz 


which  hecoiiies 


2(1111+1112) 
and  for  the  second  term 


h2(P  -mi^)  —  Acp(|  ,Ti)+Pi(mim2+3^)  — Acp(|,T]) 

o| 


(36) 


~  lim 
z— >  0 


zP^(mi+mp) 

2Jt|ii|i2 


But  this  reduces  to 


(mi+mg) 

2jtiiiH2 


P  (mi +1112)  lim 

— - z-^0 

2jtui|a2 


(37) 


and  since 


_ drii _ 

'^1  [(^2-Tlr)  (tii-^.i) 


the  second  term  in  expression  (37)  vanishes.  Finally,  therefore,  the 
vertical  induced  velocity  in  the  plane  of  the  wing  w^  becomes 
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w  =  — 
u 


2(ini+m2)y 


Acp  (|,T^)+|ii(niim2+|32)A  a9(  I,ti) 


P^(mi+iD2) 

2nni|i2 


A(p(li,T)i) 


■y*  3 


and  similarly 

for  0<]ii2<3,  0<mi<oo 


(38a) 


Wu  =  - 


2(mi+m2)</^ 


-1112“=  L 


Acp  ( |,Tl)+|i2(mini2+P^):^  Acp  (|,ti) 


St] 


P^(mi+ni2)  £ 


dll 


Acp(|i,Tii) 


27  ® 


(38b) 


and  for  P  ^102 


P^(mi+ni2)  r£ 
2rtjiiH2  JJt 


A<p(|i>ili) 


dll  dTji 


(38c) 


The  special  cases  of  equations  (38)  obtained  when  the  |,ti  coordi— 
j nates  represent  the  x^y  or  the  r,s  systems  are  given  in  the  summary 
of  this  section  (see  equations  (48)  and  (51)). 


Vertical  Induced  Velocity  in  Terms  of  the  Loading 

The  equation  for  the  loading  coefficient  in  linearized  theory  can 
be  written 


Ap  2Au 

T  "  ^ 


(39) 
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so  that  equation  (27)  for  the  vortices  yields  for  the  expression 


w 

u 


S  (iiii4-ing)zYo 
z  -^0  4jr|jLiP2 


—  +  ( ^-ii) 


1112 


(i-ni)  -i 

M-1  1^2 


+  z 


dill  (*J-0) 


The  evaluation  of  can  be  divided  into  two  steps;  first,  the  pro¬ 

cedure  necessary  for  carrying  the  derivative  through  the  first  integral, 
and  second,  the  calculation  of  I  where 


_  lim  (nii+in2)Vo 
Zr-^  0  . 

4jtpi|i2 


(ri^i)  (I-IP 


TT)2 

^2  J 


(-i-ni)  i:-  (Ml)  ^ 


+  z 


q. 


( ii^iii) 


d|^  dTli(4l) 


Again  the  order  of  integration  is  important.  To  begin  with,  tlie 
first  integration  will  be  taken  with  respect  to  rji.  Further  the  case 
0  <  mi<  p  will  be  considered.  Hence,  the  equation  for  w^  becomes 
(just  as  in  the  derivation  of  equation  (l8)) 


(mi+iite)Vo  8 

Wu  =  z-->Q  — - / 

e~>0  AnqiHg  dz  J 


^2 


mi+ma 


Lo+Li 

dTii 

Lo-Li 


(42) 


(n-^i)  —  -  (l-li)  — 


^2 


+  Z' 
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where  the  values  of  Lq  and  Li  are  given  by  equation  (17)  and  the 
limit  as  e  goes  to  zero  is  to  be  taken  first.  Equation  (36)  reduces 
to 


lim  ^Lo‘+Li‘ 

- ==  / 

e  -»0  4^111  '^Lq’-Li' 


(n-^i)  ^  + 


zm£  V 

Ell +1112 


(la^ni) 


mi+niE  J 


+  I  {^3) 


The  quantity  within  the  integral  of  the  first  term  of  equation  (43)  is 
the  same  as  the  similar  term  in  equation  (19);  hence  by  analogy 


Wu  =  - 


Vq  vP 


(i^n)  +  I 


The  evaluation  of  I  requires  some  care.  Consider  the  following 
integral  which  contains  all  the  difficulties  involved  in  I. 


T  11m  S  r  zf(yi^z)  1 


"b  >  y  >  a 


where  f(y^z)  and  its  derivative  is  boimded  and  continuous  in  the  inter¬ 
val  a  <  y  <  h. 

Integrating  by  parts 


-r  lim  S  N  .  /  X  y-1) 

-^o  “  z  -^0  T"  ^^va,zj  arc  tan - f(b,z)  arc  tan  -  + 

^  oz  z  z 


^f^(yi^z)  ^  y—yi 

- ^ -  arc  tan  — — -  dyi 

oyi  z 
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and  since 


lim 
z  »0 


^fz(yi^z) 

^yi 


arc  tan 


y-yi 

z 


dyi 


«fz(y^o)-  - 


fz(a,0)+fz(lD,0) 


this  becomes 


^  f(b)  f(a) 
y-b  y-a 


f*(yi)dyi  ^ 

y-yi 


It  is,  however,  more  convenient  to  write  Iq  in  terms  of  f(yi) 
and  not  its  derivative.  In  order  to  do  this,  it  is  necessary  to  intro¬ 
duce  another  notation  involving  singular  integrals.  The  concept  of 
Cauchy’s  principal  part  is  adopted  and  defined  in  the  following  way: 


f(yi)dyi 


y-yi 


|y-yi|  dyi 


This  procedure  can  he  generalized  (see  also  appendix  D  in  reference  4) 
so  that  for 


f(yi)dyi 

(yi-y)^ 


f(yi)d.yi 

yi-y 


l(yi)dyi 
Sy  '^a  yi-y 


(45a) 


This  definition  of  the  symbol  f  can  be  made  in  another  way.  If  the 
indefinite  integral  is  expressible  in  the  form 


f(yi)d.yi 

(yi-y)^ 


G(yi^y)  +  constant 


then 

f  =  G(b,y)  -  G(a,y)  (45b) 

(yx-y) 

Thus  the  conventional  rules  for  evaluating  a  definite  integral  from  the 
indefinite  integral  can  be  used. 

Using  equation  (45b),  one  can  derive  by  integration  by  parts  the 
relation 
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^  f(yi^O)dyi  f(b)  f(a)  f'(yi,0)dyi 


•Ja  (yi-y) 


y-b  y-a 


-I 


y-yi 


Hence  the  equation  for  becomes 


^  f(yi^O)dyi 


io=4  — ^  fz  (y.o) 

a  (yi-y) 


(45c) 


By  means  of  the  concepts  introduced  by  equations  (45),  I  can  be 
evaluated®  and  finally  Wu  can  be  written 

for  0<mi<P,0<m2<oo 

yi 


V 


2  2 

Ap  /  t  V 

—  * 


mi+iii2  /** 


4^l|l2« 


(n-^i)  ^  +  (e-ii) 


1^2 


Ap/q 


^2 


(46a) 


Similarly,  for  0<m2<P,  0<mi<co 


w„  -  2 


3i  =  _ 


mi+m2 


WiM-2 


A,./  di. 


S  -  (5-lx)  g  Ap/, 


and ,  finally,  for  p<mi  <<»,P<ni2<oo 


w. 


u 


mi4-ni2  p  p 


(^-1.)  £  Ap/q 


Vo  47tniH2j^/ 


dm 


(46b) 


(46c) 


®In  the  evaluation  of  I  the  term  representing  f(yi>z)  can  be  written 
as  a  function  of  z^.  Hence,  the  term  J(fz(y,0)  in  equation  (45c) 

vani  she  s .  _ _ _ 
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Again  the  special  cases  of  equations  (46)  vhen  "become  x,y 

or  r,s  are  given  in  the  following  summary  (see  equations  (49)  and 
(52)). 


Summary  of  Formula  for  w  and  Cp  in  the  Plane  of  the  Wing 

For  the  x^y  coordinates.—  In  these  cases  y  mi-^, 

and  m£  to . 


(x-xi)X,,(xi,yi) 


C-n  =  “  4„(x,y)  —  -  dxi/  dyi  — 

P  P  [(x-xJ^^^(y-yi)^]3/2 


(4Ta) 


=  -  ^  Ay,/ 

n  J  J  j 


dxi 


(x-xi)Xu(xi,yi) 

[  (x-xi)^-|3^(y-yi)^]®^^ 


(47b) 


n 


R  Q  2 

w^  =  -  I  Au(x,y)  -  ^  j  dxi  dyi 


A(p  (xi,yi) 


Jr  [  (x-xi)^-P^(y-yi)^]®/^ 


(4Tc) 


A'P(xi,yi) 


[( x-xi )  ( y-yi )  ^  ] 


(48) 


?  "  ”  r  — (x,y)+  ^ 


(x-xi)  (Ap/q) 


(49a) 


x-Xi)^-P^(y-yi)^ 


Wu 


dxi 


(x-Xi)  (Ap/q) 


( y-yi )  ( x-xi )  ^-P  ^  ( y-yi ) 


(49b) 


For  the  r^s  coordinates.—  .  In  these  cases  ^  y  t|  mi  —HBy 

and  ms  • 


2JtP  JtJ  [(r— ri)  (s— s,)]^/* 


(50) 
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r/Ticp  (ri,si)(iridsi 
QnJfJ  [(r— ri)(s-si) 


'  [(r-ri)+(s-si)]  (Ap/q)dridsi 
[ ( s-si )-( r-n )  ]^y  (r-ri)(s-si) 


(51) 

(52) 


It  should  be  emphasized  that  in  the  r,s  coordinate  system  the  order  of 
integration  is  immaterial. 


SOME  ALTEEKATIVE  EXERESSIONS 


It  is  sometimes  very  convenient  to  be  able  to  express  the  equations 
given  in  the  previous  sections  in  a  slightly  different  form.  Consider, 
for  example,  equation  (U)  vhich  gives  the  velocity  potential  due  to  a 
distribution  of  sources,  thus 


_  1  r  r'tfu(xi,yi)dxidyi 

where  the  area  r  can  be  defined 
( see  sketch)  as  the  area  bounded 
by  the  wing  plan  f  orm*^  and  the 
forecone  from  the  point  P(x,y, z). 
Define  now  the  area  as  being 

the  area  bounded  by  the  plane 
x=xi  and  the  wing  plan  form 
ahead  of  the  plane  x=xi  ( see 
sketch)  •  It  is  apparent  that  r^. 
is  a  pure  real  quantity  everywhere 
inside  the  area  t  and  is  a  pure 
imaginary  quantity  everywhere  out¬ 
side  T  and  inside  •  The  same 
is  true,  of  course,  of  r^®.  It 
is  clear  that  all  other  terms  in 
the  integrands  of  the  integrals 
which  have  been  considered  in  the 
preceding  section,  in  particular 
'Wu(2Ci,yi)  in  equation  (ll),  are 


The  actual  definition  of  x,  that  it  is  the  area  within  the  Mach  fore— 
cone,  is  often  replaced  by  the  one  used  here  since  the  strengths  of  the 
sources,  etc.,  are  zero  ahead  of  the  wing. 
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always  real  over  the  entire  wing  plan  forms.  Hence  the  integral 


]zj'  y^'^u(xi,yi)d.xidyi 
’’’i 


is  a  complex  quantity,  the  real  part  of  which  is  the  velocity  potential; 
thus 

1  r  r  Wu(xi,yi)dxidyi 

;  W  j  ^ 

C 

Similarly,  each  of  the  integrals  in  equations  (^7)  and  (h8)  may  he 
replaced  by  the  real  parts  of  their  values  taken  over  the  area  Tj^. 

The  evaluation  of  the  terms  involving  the  finite  part  are  peirticu- 
larly  simple  when  the  aj^ea  is  used  since,  if 


f(y)dy 

(a-y) 


=  F(y)+C 


then  for  positive  a 


and  b  ^  a 


f(y)dy 

(a-y) 3/2 


r.p.  [F(h)  —  F(o) ] 


(54) 


For  example,  consider 


y^dy 


(a2-y2)3/2 


y^dy 

(a2-y^)V2 


where  a  <  h.  From  the  relation 


y^dy 


(a2-y2) 


3/£ 


—  arc  sin 


1 

a 


together  with  equation  (54) 


I  =  r.p. 


arc  sin 


£ 

2 
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Notice  the  siM^le  extension  that 


£ 


U  (a2-y2)V2 


=  r.p, 


£ 


y^dy 


rr  —  Jf 


Turther  examples  using  the  real  part  of  the  integration  taken  over 
the  area  will  he  given  in  the  next  section  on  special  applications, 


SPECIAL  APPLICATIONS 
An  Integral  Eqviation 


Applications  of  the  results  given  by  equations  (48)  through  (52)  are 
apparent.  One  of  the  more  important  uses^  however,  comes  in  the  develop¬ 
ment  of  Integral  equations  necessary  for  the  solution  of  many  supersonic 
wing  theory  problems. 


An  example  of  such  an  application 
arises  in  the  analysis  of  the  slender 
rectangular  wing  at  an  angle  of  attack 
a.  Since  the  wing  chord  is  long  com¬ 
pared  to  its  span,  eind  since  along  the 
side  edges  the  loading  falls  to  zero, 
am  approximation  to  Ap/q  is  given  by 
the  eqxiation 

^  =  4<if(|)  (55) 

where  f(x/s)  is  an  unknown  func¬ 
tion.  The  assunqption  made  when 
usln^  equation  (55)  is,  of  course, 
to  fix  the  spanwise  variation  of 
loading  but  leave  the  chordwise 
variation  arbitrary. 


The  function  f  can  be  determined  by  the  condition  that  the  value 
of  w^  in  equation  (49a)  is  a  constant  all  along  the  center  line  of  the 
wing.  The  area  t  is  indicated  by  the  shaded  region  in  the  last  sketch 
so  that  eqmtion  (49a)  for  the  case  y=0  (and  for  added  simplicity 
3=l)  becomes  for  x  >  s 
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(f)"  -Jo" 


®  y  l-(yi/s)^d7i 


s  yi 


2  /  (x-xi)^-yi^ 


f  y  l-(yi/s)  dyi 

-( x-xi )  yiV  (2:-xi)^-yi^ 


and  for  0  <  x  <  s 


at  f  £). 


sy  Jt 


'^-<x-xi) 


Introduce  the  notation 


0,  =  e  =  |,  k.  =  ^,  k.  =  a  -  a: 


npfl  these  equations  hecome,^  since  a  =  —  ^ 


for  0  <  0  <  1 


for  1  <  0 


1  =  f(0)  +  iT  k2B2f(ei)d0i 


*  o  ^  O  P  ^ 

1  =  f(0)  +  I  /  k2B2f(0i)d0i+  I  /  Eif(0i)d0i 

O0L.1  Oq 


^The  symbols  B  and  E  indicate  elliptic  integrals.  Thus 


®n  ”, 


Bri  = 


1  /l-k/t^ 


•k  ^t^  P  ^ 

dt;  Kn  =  r  -—====>  and 
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The  solution  to  these  equations  has  been  obtained  numerically  and 
has  the  form  shown  in  the  sketch; 


Since  from  equation  (55)  the 
average-  spanwise  loading  P  can 
be  readily  calculated  as 


P 


—  dy  =  aicf 

q 


the  curve  for  average  span  loading 
can  easily  be  constructed.  This 
curve  is  also  shown  in  a  sketch 
together  with  a  portion  of  the 
variation  of  P  obtained  from  an 
exact  linearized  analysis. 

In  the  interval  where  the  compar¬ 
ison  can  be  made  (i.e. ^  near  the 
leading  edge)  the  agreement  will 
be  the  poorest  because  in  this 
region  the  spanwise  variation 
deviates  most  radically  from  the 
value  assumed  in  the  construction 
of  the  integral  equation. 


f(0) 


Drag  Eeversibility  Theorem 


The  well-known  theorem  that  the  drag  of  a  symmetrical  nonlifting 
body  is  the  same  in  forward  and  reversed  fli^t  at  the  same  speed  (see 
reference  5  or  6)  can  be  derived  in  another  way  using  the  results  of  the 
preceding  secticns. 

By  definition 

'  ^  “  s  ^  (57) 

where 

S  area  of  the  wing 
Xu  slope  of  the  upper  surface 
Cp  pressure  coefficient 


Using  equation  (47b)j 
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„  ,  ,  X  (x,y)(iy  dx 

Sjt  /  /  '‘u' 

s 


r.p 


cLki 


(z-Zi)X^(xi,yi) 


(58) 


Now  the  eqmtion  for  the  drag  coefficient  in  reversed  flow  can  be 
obtained  by: 

1.  Eeplacing  the  area  tj  by  such  that  =  S 

2.  Rotating  the  axial  system  in  the  xy  plane  through  180° 

3.  Reversing  the  signs  of  Xy^(x,y)  and  Xu(xi,yi) 

There  results 


^r  =  -s;^ 


Xu(x,y)dy  dx 


r.p. 


dxi 


(xi-x)Xu(xi,yi) 


(59) 


and  subtracting  equation  (59)  from  (58)  gives 


Cp-Cj)^  =  -  ^  r.p. J' 


xXu(xi,yi)Xu(x,y) 


r  3 


xiXu(xi,yi)Xu(x,y) 


(60) 


Since  the  symbols  xi,yi,x,y  are  dummy  variables  of  integration^  the 
last  term  in  equation  (60)  can  be  written 


dx 


_ 


and  reversing  the  operators  /dyi  /dxi  and  /dy  /dx  (but  always 
preserving  the  same  order  within  the  operation)  yields  for  the  second 
term  in  equation  (60)  the  same  expression  as  the  first  term  except  for 
sign.  Hence 


or 

Cj)  = 


as  was  to  be  shown 
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w^i(xi,yi)(ixidyi 

y  ( a+y  tan  A-xi )  ( y-yi )  ^ 

The  area  S,  being  that  of  the  wing  plan  form,  does  not  depend  on  y 
so  the  y  integration  can  he  made  first  and,  since  the  edges  of  the  wing 
are  supersonic,  the  interval  si  <  y  <  ss  must  always  contain  the  roots 
X.1  and  X,2  of  the  expression  under  the  radical.  Hence 


r.p, 


dy 


Xz 


dy 


-s  y  (3^-  tan^  A)(Xi-y)(y-j»,2)  y  (3^-  tan^  A)(Xi-y)(y-X.2) 


and  since 


then 


l 


\z  ^ 


=  It 


■1  y  (^■i-y)(y-^a) 

-4  r  r  ■w-u(xi,yi) 


y  tein^A 
An  alternative  expression  for  equation  (61)  is 


dxidyi 


(61) 


Cl  = 


y  tan^A 


(62) 


where 

tion 


a  is  the  average  angle  of  attack  of  the  surface  and  by  defini— 


a 


Vu(xi^yi) 

^o 


dxidyi 


(63) 


It  is  interesting  to  notice  that  the  lift  coefficient  for  such  a  wing  is 
the  same^as  that  for  a  two-dimensional  flat  plate  flying  at  an  angle  of 
attack  a  into  a  free  streEun,  the  speed  of  which  is  given  by  the  com¬ 
ponent  of  velocity  normal  to  the  trailing  edge  of  the  three-dimensional 
wing  just  studied.  This  result  has  been  derived  previously  in  refer¬ 
ence  7* 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  Oct.  l6,  1950. 
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